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Abstract 



We calculate the complete one-loop effective potential for SU(2) gauge 
bosons at temperature T as a function of two variables: <j), the angle as- 
sociated with a non-trivial Polyakov loop, and H, a constant background 
chromomagnetic field. These two variables are indicators for confinement and 
scale symmetry breaking, respectively. Using techniques broadly applicable 
to finite temperature field theories, we develop both low and high tempera- 
ture expansions. At low temperatures, the real part of the effective potential 
Vn indicates a rich phase structure, with a discontinuous alternation between 
confined (<fi = it) and deconfined phases (cj) = 0). The background field 
H moves slowly upward from its zero-temperature value as T increases, in 
such a way that y/gH jixT is approximately an integer. This behavior stops 
at T c = 0.722(1)^0; where fio is a zero-temperature renormalization group 
invariant scale; beyond this temperature, the deconfined phase is always pre- 
ferred. At high temperatures, where perturbation theory should be reliable as 
a consequence of asymptotic freedom, the deconfined phase (<f> = 0) is always 
preferred, and \J gH is of order g 2 (T)T. The imaginary part of the effective 
potential Vi, which originates in a tachyonic mode associated with the lowest 
Landau level, is non-zero at the global minimum of Vr for all temperatures. 
A non-perturbative magnetic screening mass of the form M m = cg 2 (T)T with 
a sufficiently large coefficient c removes this instability at high temperature, 
leading to a stable high-temperature phase with <fi = and H = 0, character- 
istic of a weakly- interacting gas of gauge particles. The value of M m obtained 
is comparable with lattice estimates. 



I. INTRODUCTION 

One of the most important features of SU(N) gauge theories at finite temperature is the 
existence of a deconfinement phase transition. The essential features of the transition have 
been well established by lattice simulations (]]]. Below the deconfinement temperature Td, 
the pressure is essentially zero, because giueball masses are large compared to T d . Above 
T d , the pressure rises, appearing to slowly approach the blackbody result for a free gas of 
gauge bosons. 
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Within the Euclidean finite temperature formalism, the deconfinement transition can 
be viewed as the spontaneous breaking of a global symmetry associated with the center 
Z(N) of the gauge group SU(N) 0. In this formalism, the temporal direction is periodic, 
with period (3 = 1/T. The Polyakov loop, defined as the Euclidean time-ordered exponential 



is the natural order parameter for the deconfinement transition. We will regard Pas an 
abstract element of the group SU (N) , and use Tr R V to denote its trace in the representation 
R. Below the deconfinement transition temperature Tj, in the confining phase, the Z(N) 
symmetry is unbroken, which in turn implies that the thermal average of the fundamental 
representation trace vanishes, (TrpV) = 0. Above Td, the Z(N) symmetry is spontaneously 
broken, and (TrpV) ^ 0. Calculation of the effective potential for V in perturbation theory 
indicates that the Z(N) symmetry is broken at high temperature, but gives no information 
about its restoration at low temperature . 

A constant chromomagnetic field is the simplest non-trivial field configuration for which 
the one-loop functional determinant can be evaluated analytically. This is essentially the 
Saviddy model |||7| . The Savvidy model is an interesting laboratory for perturbative calcu- 
lations in non-Abelian gauge theories. The essential feature of this model is the perturbative 
prediction that the vacuum of a non-Abelian gauge theory has a chromomagnetic condensate. 
As pointed out by Savvidy, this behavior can be inferred from asymptotic freedom. A direct 
path to this result is a calculation of the zero-temperature effective potential V(T = 0) 
for constant non-Abelian magnetic fields, which shows a non-trivial minimum. However, 
as first discussed by Nielsen and Olesen [Q, the zero-temperature effective potential has a 
tachyonic instability, i.e., an instability with respect to long wavelength fluctuations. This 
gives V(T — 0, H ^ 0) a negative imaginary component, indicating that a constant field 
must decay towards the true vacuum state, which is unknown. Nevertheless, the non-trivial 
minimum of V at H ^ is often regarded as evidence for the dynamic breaking of scale 
invariance in gauge theories, and for the existence of a gauge field condensate. 

The Savvidy model at finite temperature allows us to examine the coupling of the 
Polyakov loop to a gauge field condensate. Lattice results indicate that the expectation 
values of other key observables are coupled to the Polyakov loop. For example, at the de- 
confinement transition in pure SU(3) gauge theory, which is first order, both the Polyakov 
loop and the plaquette expectation values are discontinuous. In fact, the plaquette expecta- 
tion values are related to the internal energy density and pressure via the trace anomaly, and 
thus must jump at a first order transition. A simple strong coupling expansion reveals that 
plaquette expectation values depend on Polyakov loops via topologically nontrivial strong 
coupling diagrams that wind around the lattice in the Euclidean time direction. A similar 
behavior is seen for the chiral condensate @-[l]J Thus we expect to find that the field H 
couples to the Polyakov loop in the Savvidy model. 

We calculate the effective potential at all temperatures, including the effect of a non- 
trivial Polyakov loop as well as a constant non-Abelian magnetic field. From this effective 
potential, we find that the Savvidy model distinguishes between low-temperature behavior 
and high-temperature behavior in a novel manner, At low temperatures, the Saviddy model 
exhibits an intriguing oscillation between confined and deconfined phases as the tempera- 
ture is varied. We regard this as important, because understanding the mechanism which 
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determines (TrpV) as T is varied in the full theory is tantamount to understanding the 
deconfmement transition, and very likely the origin of confinement. The high-temperature 
behavior of the Savvidy model is similar to that of a free gas of gauge bosons, with a trivial 
Polyakov loop. 

Previous work on the Savvidy model at high temperatures with a trivial Polyakov loop 
has shown that the tachyonic instability found at zero temperature persists at high tempera- 
tures [12] . This is potentially much more serious than instability at zero or low temperatures, 
because asymptotic freedom is generally taken to imply the utility of perturbative calcula- 
tions in high-temperature gauge theories. The persistence of the tachyonic instability at 
high temperature is a barrier to the use of perturbation theory in this regime. Remarkably, 
a non-trivial Polyakov loop can counteract the tachyon instability at finite temperature 
|Tj|Jrj]] . We will explore this mechanism in detail in what follows. However, the basic point 
is simple: a nontrivial Polyakov loop acts as an additional positive mass term for gauge field 
fluctuations, and removes the tachyonic instability in some circumstances. Ultimately this 
approach fails: at one loop, the imaginary part of the effective potential is non-zero at the 
global minimum of V for all temperatures. An attractive possibility for high temperatures is 
that a non-perturbative magnetic mass stabilizes the model at one loop. As we will demon- 
strate below, a sufficiently large magnetic mass will not only remove the tachyonic stability, 
but also leads to a straightforward characterization of the high-temperature behavior with 
H = as the global minimum of V. 

In the process of performing the evaluation of the effective potential, we have collected 
many elementary methods for doing one loop finite temperature calculations. Many of these 
are well known. However, we have largely avoided the use of hypergeometric functions |T5|JT6 



and finite temperature zeta-function techniques [17,18), which are powerful but somewhat 



opaque. The alternative high temperature techniques we use may be of interest, independent 
of the particular features of the Savvidy model, especially since they are widely applicable, 
and retain the periodic properties associated with the Polyakov loop |TJ|. At appropriate 
points, we also stress the physical interpretation of key results, including high-T and low-T 
expansions. 

The next section develops the formalism required for the evaluation of the one loop effec- 
tive potential. Section III provides the details of the low temperature expansion, expanding 
on the results of ||13|| . A high temperature expansion is derived in section IV, and section 
V examines these results in the context of dimensional reduction. Section VI considers the 
stability of H ^ states at high temperatures. In section VII we present our conclusions. 
There are three technical appendices. 



II. THE EFFECTIVE POTENTIAL 

In this section, we introduce the formalism necessary for the evaluation of the one-loop 
effective potential for SU{2) gauge bosons at finite temperature. This is accomplished by 
calculating the partition function in background gauge, with the background field providing 
both a constant non-Abelian magnetic field and a non-trivial Polyakov loop. In order to 
carry out the calculation, we choose the color magnetic field and the Polyakov loop to be 
simultaneously diagonal. We take the color magnetic field H to point in the x% direction. 
The external vector potential can be chosen to be 
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A 2 = Hx 



T 3 



(2) 



which gives rise to a chromomagnetic field 



'12 



H 



T 3 



(3) 



The Polyakov loop is specified by a constant An field, given in the fundamental representation 
by 



Ac 



9(3 2 



(4) 



where the range of can be taken as < <fi < 2n. The trace of the Polyakov loop is given 
by 



Tr F {V) = 2cos(0/2) 
in the fundamental representation and by 

Tr A (V) = 1 + 2cos(0) 



(5) 



(6) 



in the adjoint representation. In general, the N — 1 eigenvalues of the SU(N) Polyakov loop 
are not determined by the fundamental representation trace ||20|| . However, for SU(2) and 
SU(3), the trace of the Polyakov loop in the fundamental representation does determine the 
eigenvalues. In the case of SU(2), the global Z{2) symmetry takes <ft into 2n — cf). Unless the 
Z(2) symmetry is spontaneously broken, the variable <fi must have the value it, corresponding 
to Tr F (V) = 0. 

At tree level in the loop expansion, the effective potential is given by 



V (0) = "_ H 2 

2 ' 



(7) 



the classical field energy. As explained in || , the external field H gives rise to Landau levels 
in the gluon functional determinant. The one- loop contribution to the free energy has the 



form lElO 



y{l) = V W + y(D 



V PJ (2vr) 3 



1 gH 



1 y y 1 «±L / ^£ i n 



m=0 n,± 



(3 2ix 
1 gH 



dk-> 



1 _^ 



m=0 ra.rfc 



P 2n 



dk-> 



u n - ^ ) + 2gH ( m + - ± 1 ) + k] 



P. 



1 
2 



+ 



,c„ + | ] + 2gH [ m + - ± 1 ) + k\ 



(8) 



where the uj n = 2imj(5 are the usual Matsubara frequencies, and the sum over nis over all 
integer values. The sum over m is the sum over Landau levels. The contribution Vq comes 
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from the first term, due to the neutral Z°; the second and third terms give , arising from 
the charged W + and W~ gauge bosons. The terms 2gH(m + l/2±l) are the allowed Landau 
levels for the charged gauge fields. The ±2gif term originates in the familiar coupling of 
spin to an external magnetic field. In this case, the g factor is 2, and S z = ±1. 

For m = and k$ sufficiently small, the negative sign gives rise to tachyonic modes 
which are responsible for destabilizing the original Savvidy vacuum, as first pointed out by 
Nielson and Oleson ||. As discussed in flll||14| , it is possible to avoid tachyonic contributions 
provided Trp(V) is sufficiently small in magnitude. The worst behavior occurs in the n = 
Matsubara mode and the m = Landau level. The determinant will be strictly real 
provided 

P^fgH < < 2tt - pJgH. (9) 

Thus, it is possible that the Savvidy vacuum, or a similar field configuration with apparent 
tachyonic modes, could be stabilized in the confining phase by the non-trivial Polyakov loop. 
As a practical matter, the vanishing of the imaginary part of V for this range of parameters 
provides us with an important check for various expressions. 



With the aid of the standard product representation |2~I 
cosh(x) — cos(a) 



1 



cos a 



n 



i+ 



X 



2nk + a 



(10) 



we can write the one-loop effective potential in the form: 

1 OO TT . 

Vil)=2 2 £ £ |b / 27 M cos M/^±K*0] -cos(0)} 



m=0 ± 



+ 



d 3 k 



r/3 

In {cosh [f3e k ] - 1} , 



(3 J (2tt) 3 

where Ek = \k\ and the variables uo± are defined by 

w±(m, k) = 2gH (m + - ± 1 ) + k 2 . 



(11) 



(12) 



Note that uj 2 _ can be negative. In order to obtain the correct sign for the imaginary part of 
V, it is necessary to supply the standard Feynman is prescription as needed, u 2 



l£, 



with e — > 0. The contribution of the Z can be written as 



Vr 



(1) 



d 3 k 

(2^Y 



e k + - In ( 1 



(13) 



The first term is an irrelevant vacuum energy contribution, and the second term is the free 
energy of a massless, free gauge boson. Similarly, we can write as 



vi l) 



g g_Hrdkl 1 ( 



i.n,l 



-f3(jJ±—i(: 



(14) 
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There is a clear separation of the zero-temperature part of the effective potential and 
the finite temperature part. We write the finite temperature part of V and V± as Uq and 
U±, respectively. All ultraviolet divergences reside in the zero temperature part 

V<V(T = 0) = V {1) (T = 0) + V2\T = 0) 



m=0 ± 



The finite, if-dependent part of V^ 1 -* (T = 0) was first derived by Nielson and Oleson || ; for 
the reader's convenience, we review their derivation in Appendix 1. With an appropriate 
choice of renormalization constant, can be absorbed into V^'(T = 0) and V(T = 0) is 
given by 

V(T = 0) = V {0) + V (1) (T = 0) 

llg 2 H 2 , /VfA (off) 2 , , 

where fio is a renormalization group- invariant parameter that sets the scale for the gauge the- 
ory. The minimum of the effective potential occurs at gH = A|, where As = /io ex P ( — 1/4). 
Note that the combination gH is renormalization group invariant in background field gauge, 
and henceforth g and H will generally appear together. 



III. LOW TEMPERATURE EXPANSION 



The one-loop finite temperature contribution to V from the Z° gauge boson, which we 
write as Uq, is given by 



Un 



2 



d 3 k 
(2^ 



In ( 1 



-2- 



7T 



90p 4 



(17) 



This is the free energy density of a boson gas with two spin degrees of freedom. The phase <j) 
does not appear because the Z° is charge-neutral. Because the free energy is extensive in the 
volume, this is also the negative of the Z° contribution to the pressure. The W ± contribution 
U± is somewhat more difficult to evaluate. Expanding the logarithm, the contribution of 
the Landau levels to the effective potential can be written in the form: 



2 " gH r dk 
~?h± 2vr J 2n 



cosine 



-n(3uj±(m,k) 



n=l 



n 



(18) 



This expression has a natural interpretation in terms of path integrals. At finite temperature, 
there are particle trajectories which wind around space-time in the Euclidean temporal 
direction, and are thus topologically non-trivial. For such trajectories, there is a factor of 
exp (in<fi) when the net winding number is n. Alternatively, we can consider i(p as a chemical 
potential continued to imaginary values. 

Some care must be exercised, because uj-(Q,k) = ^2gH ^0 + ^ — 1 
for k sufficiently small. The contribution to U± from u_ with m = is 



A; 2 is imaginary 
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2 gH ~ cos(n0) /■ _ n 



i/3tj_(0,fe) 



2tt 



(19) 



We define k 2 = gH so that u>~(0, k c ) = 0. Then u(—, 0) can be written as 



u(- 0) 



2^i7 y, cos(ra<; 



o 2^ 



e -n/3a._(0,fc) + 



00 d& 

fc c 27T 



-n/3a)_(0,fe) 



(20) 



The analytic continuation of the first integral is carried out using the general prescription 
that —is is added to 2gH (n + | ± 1 j , analogous to the standard m 2 — is prescription. For 
the case of the (— , 0) contribution, this means that gH —>■ gH + is. The first integral can 
be written as 



o 2vr 



tt/2 

— cos 9 [cos (n(3k c cos6)} + — -iriJi(nf3k c ) ■ (21) 

27T 47T 



while the second integral is |22 

^ dk 



I" — e-W* 511 *! = -^(n^jfcc) - — T /2 ^ cos (n/3ife c cos 0) cos 0. 

Jfc c 27T 4 27T JO 

When the two terms are added, the remaining integrals cancel, yielding the result 

(gH) 3/2 00 



(22) 



u(- 0) 



E 



cosine 



27T/3 ^ n 
The other terms all have the form 

2 gH ™ cos(n0) r dk 



Y 1 (nf3k c )-iJ 1 (nf3JgH) 



(23) 



2 gH - cos n0 y 



2^ 6XP 



where m > lfor u(—,m). Using the integral representation 



-n(3\j2gH (m+-±l) +P 



(24) 



Ki(n(3M) 



2M 



dk e 



-n/3Vk 2 +M 2 



(25) 



we obtain 



it(±, m) = — 



gH_ 



l2gH[m+-±l) £ 



cos m 



n=l 



77 



#1 



n/3 1 /2^(m + -±l 



(26) 



The contribution of each w(+, m) is doubled by the u{— , m+ 1) term, except w(— , 1) stands 
alone, giving a contribution to U± of the form 



3/2 



E 



cosmc 



n 2 P ^ n 



Kx(npJgH) + 2J2 V2m + 3K 1 (npJ(2m + 3) £#) 



m=0 



(27) 



Combining the results for V(T = 0), Uq, and U±, we obtain finally a renormalized effective 
potential with real component 
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v R 



n(gH) 2 (gjT 
48tt* n 



2tt 2 
90^ 



(gH) 



3/2 oo 



E 

n=l 



3/2 oo 



cos nc 



7T 2 /3 

and imaginary component 



E 

n=l 



n 



cos nc 



K 1 (n/3y/g~ff)-^Y 1 (nPy/g~H) 



n 



J2 V2^T3K 1 [np^J(2m + 3)gH] 



(28) 



m=0 



V/ = - 



(^) 3/2 ^ cos(n0) jf-r-jn 



(29) 



n=l 



Note that Vr = when H = and T = 0. Numerical testing verifies that V/ is 
indeed zero whenever (3y/gH < <p < 2tt — fiy/gH as required from the expression in terms 
of Matsubara frequencies. A similar calculation of V^and V/ has also been performed by 
Starinets, Vshivtsev, and Zhukovskii []T4j| . They also noted the condition for stability, Eq. @. 
In their work, it appears that the Bessel functions J\ and Y\ were inadvertently interchanged 
in the formulae for the real and imaginary part of the potential, but otherwise our formulae 
are identical. Our results are in numerical agreement with the exact result for Vj derived by 
Cabo, Kalashnikov, and Shabad [23] for the case = 0, and with a generalization of their 
expression for ^ 0, derived below. 

At low temperatures such that f3 2 gH -C 1, the dominant contribution to Vr comes from 
the terms involving Y\, which arise from the tachyonic mode, and are entirely responsible 
for the temperature-dependent part of Vj. Thus, we may write at low temperatures 



>9H 



2tt Jo 
The real part is 

ReU± w 



dk 
2^ 



gHl 



(30) 



/gH dk 
2vr (3 Jo 2^ 



In 2-2cos l/3JgH-k 2 + 



+ ln(2-2cos [PyfgH 



k 2 



(31) 



and the deriviative of ReU± with respect to is 







ReU± 



gH 1 fvP dk 
~2^p Jo 2^ 



sin 



gH - k 2 



cos 



gH -k 2 + 
sin (j3y/gH - k 2 - . 
- cos (Py/gH - k 2 ■ 



(32) 



It immediately follows that the minimum of U± must occur near either = or = it 
for low temperatures. Careful numerical analysis of the low temperature formulae confirms 
that the global minimum indeed is given by either = or = tt. 
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A similar analysis of the derivative with respect to gH shows that the minima of ReU± are 
found at (3y/gH ~ 2n7r when = 0, and (3^/gH « (2n + 1) tt for = ir. Numerical analysis 
shows that = is the global minimum of Vr for all temperatures T > 0.722(1)/xq- As the 
temperarature is lowered, the global minimum of Vr alternates between = and = tt. 
For the range of temperatures 0.198(l)/xo < T < 0.722(1)//q, the global minimum is at = 7r, 
(3^/gH m tt. As the temperature is lowered towards zero, the minimum of gH moves toward 
its zero-temperature value, and the global minimum of Vr continues to alternate between 
= and = 7r, with corresponding changes in gH. In figures 1 through 4 we plot Vr as 
a function of \J gH jixT for = and = n for successively lower temperatures, showing 
the alternation of the minima. Note how the minima occur at integer values of the variable 
y/gH/itT. Figure 5 shows how a rapid oscillation is superimposed on the zero-temperature, 
0-independent behavior of the Vr at very low temperatures. Numerical computation of Vi 
shows that Vj at the global minimum for any given temperature. Thus the Savvidy 
model is unstable at low temperatures. However, it is striking that the low-temperature 
behavior of the model shows a strong role for the Polyakov loop, with many local minima 
ofVR. 



IV. HIGH TEMPERATURE EXPANSION 

In order to develop a suitable high temperature expansion for V, we return to the original 
expression for v£\ using the standard device of Schwinger to express the logarithms of 
equation (§) as an integral 



-EE 



1 gH r dk 



m=0 i 



^ P 2tt J 2vr Jo t 



dt 



exp 



2im — 0\ / 1 \ i 

J +2gH(m+-±l) k 2 



where we have used the symmetry of the expressions under n < — 
contributions from W + and W~ . We apply a 6* 4 identity of the form 



E ex P 



j 2 (<t>-2nnf 
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E ex P 



gg 2 
At 



+ ip4> 



(33) 

nto combine the 
(34) 



and perform the trivial integration over k to obtain 



exp 



gp 2 
At 



+ ip4> 



(35) 



In the zero temperature limit, only the 0-independent p = term survives. We henceforth 
omit this term, since it is included in V^'(T = 0), evaluated in Appendix A. Using the 

symmetry under p < > —p, the remainder of the above expression is the finite temperature 

contribution 



m=0 ± p=l ^ H JU 1 



-t ■ 2gH [m+ - ± 1 



exp 



At 



cos(p0). (36) 



We perform the sum over m and ±, obtaining 
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qH f°° dt tnH ~ 
4tt 2 Jo t 2 



P 2 p 2 
At 



COs(j90) 



(37) 



where the second term is due to the unstable mode. This integral must be defined by analytic 
continuation and is responsible for the imaginary part of U±. After the shifts t — > t/gHior 
the first term and t — > —t/gHior the second, we have 



U± = 



(gHf r^dt 

An 2 Jo t 2 



e 1 coth(t)J ex P 
P =i 



(3 2 gHp 2 



+ n / — e *]Texp 
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The first term can be decomposed using the expansion 

< 00 
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where B n is the n'th Bernoulli number. We define 
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The unstable mode gives a contribution to U± which we write as 



J-gHf r dt_™ 

f - 1 = + -^Jo ¥ e g exp 



At 



COs(p0) 



(42) 



(43) 



using the analytic continuation discussed in section III. The complete expression for U± is 
the sum 



U± = f + fl + fk>2 + f-L 

Using the integral representation of the Bessel function 
^ = 2 I2J Jo dt ¥^^ 



z 

'-At 



(44) 



(45) 
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we write fo as 



2gH " cos(p0) r— 



(46) 



and /x is similarly 



/i 



Git 2 



(47) 



P =i 



The /_i term can be done by considering the integral as a function of —gH, and gives 

2tf(-Stf) 1/2 ^cos(W>) 



7T 2 /3 



E 

p=l 



P 



-K^ppJ-gH) 



(48) 



where the analytic continuation is again specified by lim e ^ V~gH — is = —i^gH. These 
intermediate forms will require resummation to obtain the high temperature limit. 

Resummation of f and f\ is made possible by a set of identities for Bessel function 
sums which are quite useful in finite temperature field theory [Jj"9"| . For completeness, a brief 
derivation is given in Appendix B. The first of these identitites is [ETf] 



J2 K o(pz) cos 
p=i 



x '7T 

+ ? 2 



z 2 + (0-2Ztt) 2 27r Kl 



(49) 



where 7 is Euler's constant. The notation J2i indicates that the singular l/|Z|term is 
omitted when I = 0. This leads immediately to the formula 
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The remaining Bessel function identities yield 
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and 
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(52) 



Although both / and /_i appear to contain polynomials in which are not manifestly 
periodic, these terms are the representation on the range Oto 27rof periodic functions. As 
explained in Appendix B, they are obtained from the Bernoulli polynomials. 
The analytic continuation of the logarithm in /_i gives 
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Note that imaginary terms can potentially arise from a finite number of the square roots in 
depending on the value of 0. 
The remaining term, fk>2, is 
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Once more performing a # 4 transformation we obtain 
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The first term in curly brackets represents a finite contribution to the renormalization of 
{gH) 2 and is evaluated in Appendix C; the integral over t in the second term can be per- 
formed analytically, and we obtain 
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The numerical value of the constant is approximately C\ = —0.01646. 
The complete expression for U± is 
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This can be added to the previous, much simpler, expressions for V(T = 0) and U to give 
a complete expression for V. With the exception of the first term, the entire expression is 
manifestly periodic in 0. This first term is essentially the fourth Bernoulli polynomial, and is 
valid as written for the range < < 2n. Note that the order T 2 terms involving the second 
Bernoulli polynomial have disappeared from the final expression, the result of a cancellation 
of contributions from f and The omission of the tachyonic mode contribution /_i led 

0)0 



The 0(T 4 ) term dominates 
) is always zero at the global 



to a spurious T 2 term in an early calculation of V(<f> 
the effective potential at high temperatures, which implies that 
minimum of the effective potential. 

It is possible to extract from V a simple representation for Vj. Let L + be the largest 
integer such that 2nL + < /3y 'gH + and L_ be the smallest integer such that 27rL_ > 
—/3y/gH + <p. Then Vj is given by 
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which represents a generalization of an expression first derived by Cabo et al. p3| for the 
case = using different methods. In addition to the derivation from the high temperature 
expanion, we have also derived this result for arbitrary <fi using their methods, and have 
verified numerically that this expression is equal to the low-temperature form derived in the 



previous section; see reference [13| for graphs of this function. 

As a simple check on our results, we consider the H — > limit. This limit follows quickly 
from the expressions for V(T = 0), Uq, and U±. We have 
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and of course Vi(H = 0) = 0, in agreement with the results of references It is invariant 

under the substitution <p — * 2tt — 0, reflecting the Z(2) invariance of the gauge theory. The 
minimum of Vr occurs at <fi = 0, or equivalently <fi = 2tt, where Trp(V) = 2. For <p = 0, Vr 
is simply the free energy of a black body with 6 degrees of freedom, resulting from 3 colors 
in the adjoint representation, each having 2 spin states. This is the naive behavior expected 
at high temperatures: a free gas of gauge bosons. 

The = limit is somewhat more complicated than the H = limit. The complete 
expression is 
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Let L be the largest positive integer such that 2nL < (3\fgH. Then the imaginary part of 
V, Vi , is given by 
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which is precisely the result of Cabo et al. [2vJ. The high temperature limit of Vr to order 
T°is 



-6- 



7T 



11 



24tt 2 
W /2 

27T/3 



y 
In 



90/3 4 
4tt i 



+ 7 



1 ^ 2 2fe 5 2fe 



6 + 27rV2 E 



fc=2 



(2*)! 



32vr 2 



T{2k - 3/2) 



+ 0((3 2 (gHf 



(62) 



The sum in the last term is converted to an integral and evaluated numerically in Ap- 
pendix C. The result is 
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where C 2 has the approximate value C 2 = 0.82778, in agreement with the work of Ninomiya 
and Sakai |12j and of Persson [lffi|. Appendix C also proves that our expression for C 2 is 
equivalent to that given in reference [12]. In the high temperature limit, only the L = 
term contributes to Vj, which gives 
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which also agrees with references [12.25 
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We can identify the coefficient of {gH) 2 

ah CO = 



in V R as l/2g 2 eff (T), 



i.e., 
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which is positive because /3/xo < 1 and goes to zero as T 
freedom. The minimum of Vr occurs at 
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The minimum value of this dimensionless variable goes slowly to zero at T goes to infinity. 
As T ranges from 2/zoto 10/io, the minimum decreases from 0.248 to 0.131, so the assumption 
that f3 2 gH is small can be justified at temperatures not much larger than /io. Thus perturba- 
tion theory requires that the Saviddy model has H ^ at arbitrarily high temperatures, and 
the standard perturbative state (H = 0) is a local maximum of the effective potential. This 
result must be considered more reliable than the similar zero-temperature result, because 
asymptotic freedom applies. However, the free energy of the H ^ state continues to have 
a non-zero imaginary part at high T, so the Saviddy state (constant H ^ 0) is also unstable. 



Our results here are in complete agreement with the earlier work of Ninomiya and Sakai |12 
The deviation of V^trom the black body, (H = 0) result is less than 1% for all temperatures 
above 2/zo, so any indication of this effect, if it exists, is essentially unobservable in lattice 
determinations of the pressure and other thermodynamic quantities. 



V. HIGH T BEHAVIOR AND DIMENSIONAL REDUCTION 



As shown in the previous section, the leading behavior of V at high temperature is 
of order T 4 , and is independent of H. Due to the cancellation of order T 2 terms, the 
next-to-leading term is of order T. The origin of this 0(T) term has been obscure, and 
was non-trivial to obtain via zeta function methods even in the case of free fields [ I7LIT5 . 
In this section, we show how the 0(T) term arises naturally from the n = mode in 
the context of dimensional reduction fl26[-gj|. The application of dimensional reduction is 
straightforward: the functional determinant can be regarded as an infinite product of three- 
dimensional functional integrals in which each Matsubara frequency n has a mass of \2im/(3\. 
The n^0 modes contribute a T 4 term which is independent of H, as well as a logarithmic 
correction term to the classical action. The n = contribution must be 0(T) because the 
only T dependence in this mode arises from the replacement 
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The contribution of the n = mode to V, which we write as V n= o , can be treated by the 
techniques developed above. After using Schwinger's proper time representation, V n =o has 
the form 
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which becomes, after summation over m and integration over k^, 
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which may be compared with the similar terms contained in Eq. [)?]. The loss of periodicity 
in is expected when treating only the n = mode; similar behavior has been observed 
in calculations of the dimensionally reduced theory with H = [30fl. Note that the Z(2) 
symmetry does remain as the discrete — > —0 symmetry. 

Direct comparison with the high temperature expansion for V is simplest when 0is 
taken to be 0, which is appropriate at high temperatures. A different form for V n= o can be 
obtained in this case by integrating first over k% and then over t, deferring the summation 
over m. This gives 
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where the Hurwitz zeta function is defined via the series 
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The real part of this expression has the same form obtained by Persson . The approximate 
value, after applying a reflection formula for ( (—1/2, 1/2) and summing the resulting series 
numerically, is 
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This result for the 0(T) part of V is in exact agreement with other results in the case = 0, 
but here is completely attributable to the n = Matsubara mode. The full mode sum is 
necessary to recover periodicity in 0. 



VI. MAGNETIC MASS AND STABILITY OF H ^ AT HIGH TEMPERATURE 



It has been known for some time [12|] that the Savvidy model at high temperature 
continues to exhibit the pathologies associated with its zero-temperature behavior. The real 
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part of the one- loop effective potential has a minimum at H ^ 0, and the imaginary part 
of the potential is non-zero at that minimum. This result must be considered more reliable 
than the similar zero-temperature result, however, because asymptotic freedom applies. As 
we will show below in the context of the Savvidy model, the effect of any gluon condensate 
at high temperature would be difficult to observe in lattice calculations of thermodynamic 
quantities, because the condensate effects are of order g 6 T 4 . Nevertheless, it would be 
troubling if finite temperature effects did not act to eliminate the gluon condensate at high 
temperatures. It is intuitively appealing that the high temperature behavior of a gluon gas 
approaches that of a non-interacting relativistic gas. If the high temperature behavior is 
fundamentally this simple, the Sawidy instability must be removed by some mechanism. 

The resolution of the high-temperature stability issue is also of interest on phenomeno- 
logical grounds. It has been suggested by Enqvist and Olesen [|3T| that the large-scale 
non-Abelian magnetic fields in the early universe may have provided a mechanism for seed- 
ing the galactic dynamo. In their work, they used an approximation to the Saviddy model at 
high temperature. However, certain of their assumptions were questionable. In particular, 
they assumed that gH would remain near its T = value at high temperature. As we have 
seen, the log (H) term responsible for H ^ at zero temperature is replaced by a log (T) 
term at high temperature, and the value of H at high temperature need not be commensu- 
rate with the zero temperature value. They also assumed that the electric screening mass 
M e , which is of order gT, plays a role in overcoming the tachyonic instability. As we discuss 
below, it is rather the magnetic screening mass M m , believed to be of order g 2 T, which is 
relevant. Later work by Elmfors and Persson |32|] used the magnetic mass rather than the 



electric mass. However, they also assumed that a spontaneously generated magnetic field 
at high temperature would be close to the zero temperature value. Using a renormalization 
group argument, they showed that gH at T = is always less than M^, and concluded 
that a spontaneous magnetic field would be irrelevant at high temperatures. They therefore 
considered only the case of an externally imposed field in detail. 

An ambitious analysis of the behavior of H at high temperatures has been attempted in 



the recent work of Skalozub and Bordag |33| , which includes the effect of two-loop and ring 
diagrams. They assume the magnetic screening mass originates solely from having H ^ 0. 
Their result for the real part of the effective potential indicates that a non-zero H is favored 
at high temperature, with gH of order g 8 ^ 3 T 2 . However, they also note that their expression 
for the imaginary part of the effective potential is non-zero, again implying that H ^ is 
unstable. Unfortunately, our results for the free energy at one loop disagree with theirs in 
the term proportional to T. There are additional order T terms in our expressions for f\ 
and fk>2 that are responsible for the difference. Our results are in agreement with the earlier 



calculations of Ninomiya and Sakai |L2| and of Persson p5 |. 

Rather than exploring a specific origin for the magnetic screening mass, we simply assume 
that a magnetic screening mass M m of order g 2 T is present, and consider the consequences 
of that assumption. As we will show below, if the constant of proportionality is sufficiently 
large, the Savvidy instability is removed, and H = is favored at high temperature. The 
addition of the magnetic mass requires the replacement 

(74) 
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in all sums over Matubara modes. We need only consider the impact of M, m on the n = 
Matsubara mode, because is assumed of order g 4 T 2 , and is therefore negligible compare 
to the T 2 term occuring when ti^O. The contribution of the n = mode to the effective 
potential is 
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appropriate at high temperature, we have in the case < gH 
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0, as is 
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and 



V, 



n=0 



(gH) 



3/2 r 



E 



2 2k B 2k T(2k - 3/2)({3 2 gH) 2k ~ 3 / 2 



4^/2/3 (2k)\ ((3 2 M 2 m + (3 2 gH) 



2fc-3/2 



+ 



gH_ 

2np 



M 2 



gH (77) 



for the case > gH. Using the techniques developed in Appendix C, the evaluation of 
the infinite series can be transformed into the evaluation of an integral. We find numerically 
that V n= o can be well-approximated for all values of the magnetic mass by the k — term 
in the sum. This approximation is worst at M m = 0, for which the error is still less than 
20% ; the error falls to about 6% for = gH, and approaches zero for ^> gH. Thus, 
for M^j < gH, we write 
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with a similar expression in the case > gH. All other contributions to V are as before. 
As in our previous, perturbative analysis of the high-temperature behavior of V, we can 
drop all terms in the dimensionless potential (3 a Vr which are of order ([3 2 gH) 3 or higher. 
Then (3 a Vr can be written as 
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where we have again used the effective, temperature dependent, coupling constant defined 
by 
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We now assume that M m can be written to leading order as M„ 
x — (P 2 gH) / giff- We then have 
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Since only the term in square brackets need be minimized, this form explicitly shows that the 
minimum value of g H will be of order gjffT 2 , and that the free energy consists of the usual 
T 4 blackbody term plus a gtffT A correction due to chromomagnetic effects. Of course, this 
expression does not include those perturbative higher-loop corrections to the free energy, 
which start at order g 2 f fT A but are independent of H El, nor does it include higher- loop 



corrections depending on H, which start at g e $fT [[35 
For c > l/vrv^ ~ 0.225, the non-trivial minimum of 



-x 2 - — (c 2 + x) 3/2 + 6(c 2 -x) — Vc^ 
2 3vr V / V / 2tt 



X 



occurs at x = c 2 . This is degenerate with x = when 



Ccrit 



2 

37 



(2) 



3/2 



1 



0.388. 



(82) 



(83) 



If cis larger than c cr u, then H = is the minimum of Vr and Vj vanishes. In other words, 
the one-loop prediction is that H = provided 
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The numerical coefficient is likely to be changed by higher-loop effects, but indicates that a 
sufficiently large magnetic mass will lead to H = 0. 

In principle, the magnetic screening mass can be determined from lattice measurements 
of the gluon propagator. In Landau gauge, Heller et al. found that M m was well fit by 
0.456(6) g 2 (T)T over a wide range of temperatures [36]. In maximal Abelian gauge, Cuc- 
chieri et al. determined M m = 1.48(17) T at T = 2T C , which is equivalent to 0.505 g 2 (T) T 
using the one-loop form for the running coupling constant assumed by the authors J37| . 
Later, more extensive, work by the same authors [3~8~| found a complicated gauge- and 
volume-dependent structure in the magnetic propagator at low momentum inconsistent with 
a simple pole mass. Further progress in extracting a magnetic gluon mass from lattice simu- 
lations is thus dependent on progress in understanding the low-momentum structure of the 
finite-temperature gluon propagator. Although values for M m obtained from lattice simula- 
tion are large enough to possibly argue against a Savviddy instability at high temperature, 
the relatively small difference, i.e., c CT u = 0.388 versus 0.456, combined with many theo- 
retical uncertainties, provide no definitive resolution of the stability issue. Indeed, the fact 
that the values are commensurate may indicate that the generation of the magnetic mass is 
intrinsically related to the Saviddy instability in some unknown way. 



VII. CONCLUSIONS 

The ground state of non- Abelian gauge theories has two related and important features: 
confinement and scale symmetry breaking. The Savvidy state, consisting of quantum fluc- 
tuations around a constant chromomagnetic field at zero temperature, provides some insight 
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into the nature of scale symmetry breaking. As we have seen, the analysis of quantum fluc- 
tuations around a constant field at finite temperature in an SU(2) gauge theory allows us 
to study aspects of the interplay between scale symmetry breaking and confinement. 

At low temperatures, a complicated behavior emerges from the one-loop effective poten- 
tial. The real part of V, Vr has minima near integer values n of the dimensionless variable 
^J[gH)/iiT. As T goes to 0, the global minimum corresponds to higher values of n, and H 
approaches its zero temperature value. The preferred value of alternates discontinuously 
between and tt, representing an increasingly rapid set of transitions between confined and 
deconfined phases as the temperature approaches zero. However, the imaginary part of the 
effective potential, Vj never vanishes at the global minimum of Vr. Therefore the Savvidy 
model is unstable at low temperatures. While there is no reason to trust a one loop pertur- 
bative calculation at low temperatures, these results demonstrate that gluon propagation in 
a non-trivial background can lead to confinement at low temperatures. 

At sufficiently high temperatures, the leading term in the free energy, which is propor- 
tional to T 4 , demands that = 0. The leading behavior of Vr is that of a free gas of 
massless, non-interacting gluons. The dominant subleading contribution to Vr at one loop 
comes from the n = Matsubara frequency. This term leads to gH being of order g 4 T 2 , 
in turn making a contribution of order g 6 T A to the free energy. The imaginary part of V 
remains non-vanishing, and thus a constant chromomagnetic field is unstable in perturba- 
tion theory at high temperatures. However, a magnetic mass of the form M m = cg 2 T would 
naturally alter the one-loop results in such a way that <ft — 0, H = is favored in the high 
temperature limit, provided that the coefficient c is sufficiently large. This would restore 
the picture of the high temperature state as a plasma of weakly interacting gluons. In this 
case, the magnetic sector would still contribute to V at order g 6 T 4 . 

The existence of a sufficiently large magnetic mass is consistent with our knowledge of 
the magnetic sector gleaned from the work of Karabali and Nair on three-dimensional gauge 
theories in the Hamiltonian formalism [|39| -|42"|1 ■ They are able to extract a three-dimensional 
gluon mass, corresponding to a magnetic mass in four dimensions. Their estimates of this 
mass are consistent with those obtained from lattice simulation. As we have seen, the masses 
obtained in the case of SU(2) are not sufficiently large to confidently assert that H = is 
stable. It would be very useful to incorporate the results of Nair et al. into an effective 
potential for the finite temperature, four-dimensional theory. 

Over the years, many kinds of field configurations have been suggested as being re- 
sponsible for confinement. A constant chromomagnetic field is the simplest non-trivial field 
configuration for which the associated functional determinant can be obtained analytically. 
It is interesting to speculate as to what features of this model might carry over to more 
complicated field configurations, and how the Polyakov loop might be driven to confining 
behavior. The Savvidy state does appear to distinguish at one loop between the low tem- 
perature regime, where (TrpV) is a sensitive function of the temperature, and the high 
temperature regime, where (TrpV) goes to its maximal values of ±2 in SU{2). Perhaps this 
behavior is some distant relative of the deconfinement transition. 

Recent work on phenomenological models of the deconfinement transition may provide 
some direction for further theoretical investigation. In collaboration with Travis Miller, 
we have constructed two models for the free energy which give rise to confinement at low 
temperatures f2(J. These models confine by adding a phenomenological non-perturbative 
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term to the free energy which depends on the Polyakov loop eigenvalues. Both models 
account very well for many features of the deconfinement transition observed in lattice 
simulations. A similar approach has also been developed by Dumitru and Pisarski |H3-E5|. 
Based on their work, Sannino has recently proposed a phenomenological effective potential 
for finite temperature QCD in which the Polyakov loop is coupled to the scalar glueball 
field [j46| . Because the operator Tr couples to the scalar glueball sector, our expression 
for the Sawidy effective potential shares certain features with this more phenomenological 
approach when H 2 is identified as the glueball field expectation value. In either approach, 
however, it does not appear that the glueball field is driving the deconfinement transition; 
some other mechanism is still required. 

A complete field-theoretic description of confinement, perhaps using other observables in 
a role similar to H in the Sawidy model, should produce an effective action for the Polyakov 
loop which yields confinement at low temperature in a natural way. 



APPENDIX A: ZERO TEMPERATURE 



In this appendix, we derive the zero temperature part of the effective potential, following 
the approach of reference |§. The only if-dependent part is V$p(T = 0) which we can write 

as 
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after analytic continuation of the k integral to 1 — e dimensions. We have introduced an 
i8 to correctly treat the imaginary part. This introduces the scale parameter fi. Performing 
the k integration, we obtain 
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where the generalized Riemann function ((s, a), or Hurwitz function is defined by 

OO 

The term in square brackets can be written as 
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We may evaluate ((— 1 + |) to order e using formulas from Actor |L7|,[T8|. To order e 
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where C is a real, finite constant. Using 
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we find 
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Renormalization removes the 1/e pole and the renormalization coupling g(fi) is defined in 
the standard way. The complete zero-temperature effective potential is 



V(T = 0) = V {0) + V (1) (T = 0) 
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We are free to introduce a new renormalization group invariant /j,q such that 

I2tt 2 
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yielding finally 
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APPENDIX B: IDENTITIES FOR BESSEL FUNCTION SUMS 



In this appendix, we prove the Bessel function identities: 
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The right-hand sides of the second two equations are valid for <ft in the range < <p < 2tt. 
They can be extended to all real values if 0is replaced by \<p\ mod 2tc on the right hand side 
of the equation. The second and third identitites follow from the first 
Using a standard integral representation |H] 
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00 cos(pt) 



lo \/t 2 + z 2 ' 
we can write the left-hand side of the first identity as 
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We introduce a regulating mass /i, which will be taken to zero at the end of the calculation, 
and add and subtract the divergent m = term, obtaining 
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Using the Poisson summation technique in the form J2m $ i^z — m) = J2n ex P {—2irink z ), we 
obtain 
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The first integral gives a sum of screened Coulomb, or Yukawa, potentials 
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and both terms appear to be problematic as \i — > 0. The first term can be made finite in 
this limit by subtracting the contribution at r = = for ti ^ 0. Using the notation yj n 
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to denote a summation over all n with the omission of the singular term when n = 0, we 
have 
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The first sum over n is finite as /x — > 0, The second sum over n can be explicitly summed 
and the final integral carried out, leading to 
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In the limit /i — > 0, we finally obtain 
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The second identity is obtained from the first using the identity 
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, -K„(z) = -K^(z) - -K„(z). 
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It follows immediately that 
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so that 
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p=i f p= 
where C ((f)) is an unknown function to be determined later. Integration yields 
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The function C(<f) is determined using the behavior of K v {z) for z — > 



and the standard result 
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valid for < <fi < 2tt. It can be extended to all real values if <j)is replacd by \(f>\ mod 2tt on 
the right hand side of the equation. This implies the leading behavior of the sum as z — > is 
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giving us finally 
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The third identity is obtained from the second by a similar argument, combined with 
the identity 
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also valid for < 6 < 2it. 



APPENDIX C: SUMS INVOLVING BERNOULLI NUMBERS 



In this appendix, we show how to evaluate the constants C\ and C 2 introduced in section 
IV. Using the generating function for the Bernoulli numbers 
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Using the same method, we can also evaluate 
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Ninomiya and Sakai give an alternative expression for this constant 
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the individual terms of which are very badly behaved. Writing this expression 
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one can show that the difference is zero using integration by parts. 
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FIGURES 
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FIG. 2. Vr versus the scaled variable yfgH/nT for T = 0.15/x 
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